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Abstract, he quantum cat map is a model for a quantum sys- 
tem with underlying chaotic dynamics. In this paper we study the 
matrix elements of smooth observables in this model, when tak- 
ing arithmetic symmetries into account. We give explicit formulas 
for the matrix elements as certain exponential sums. With these 
formulas we can show that there are sequences of eigenfunctions 
for which the matrix elements decay significantly slower then was 
previously expected. We also prove a limiting distribution for the 
fluctuation of the normalized matrix elements around their aver- 
age. 



1. Introduction 

The quantum cat map is a model for a quantum system with under- 
lying chaotic dynamics that was originally introduced by the physicists 
Hannay and Berry [TU]. This model can be used to study the semiclas- 
sical properties of such systems [21 [6J EEl HE] . The classical dynamics 
underlying this model is the discrete time iteration of a hyperbolic map, 
A G SL(2,Z), on the torus, T 2 = R 2 /Z 2 . In order to quantize the cat 
map, for every integer N (playing the role of the inverse of Planck's 
constant) the Hilbert space of states is T-Ln = L 2 (Z/iVZ). For every 
smooth real valued function / there is a quantum observable, i.e., a 
Hermitian operator Op JV (/) : Hn — ► T~Cn- The quantum evolution is 
given by a unitary operator Un(A) on Hat. 

For generic quantum systems with underlying chaotic dynamics, it is 
believed that matrix elements of smooth observables tend to the phase 
space average of the observable in the semi classical limit. In order to 
test this phenomenon in the quantum cat map model, Kurlberg and 
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Rudnick introduced hidden symmetries of this model, a group of com- 
muting operators that commute with Un(A), they called Hecke oper- 
ators |15] . They showed that for any sequence of Hecke eigenfunctions 
(i.e., joint eigenfunctions of all Hecke operators), the corresponding 
matrix elements converge to the phase space average as iV — > oo. To 
be more precise they showed [T5J Theorem 1] that for any / e C°°(T 2 ) 
and ip G T~Cn a Hecke eigenfunction the matrix elements satisfy 

|(o Pjv (/)v,V)- / /| </, e A^ +e . 

They remarked [TSJ Remark 1.2] that the exponent of | is not optimal 

and that the correct bound should be 0(N~^ +e ), in accordance to the 
second and fourth moments. For N prime (and consequently also for 
N square free) this is indeed the correct bound [U [9]. 

Remark 1.1. We note that without the arithmetic symmetries these 
bounds hold only if the spectral degeneracies are sufficiently small. In 
fact, there are sequences of eigenfunctions (where the degeneracies are 
exceptionally large) that don't converge to the phase space average at 
all. For these eigenfunctions the matrix elements localize around short 
periodic orbits in the sense that the coresponding limiting measure 
contains a component that is supported on the periodic orbit [6]. 

In [TB] Kurlberg and Rudnick went on to investigate the fluctuation 
of the normalized matrix elements, 

(1.1) Fj N) = ((OvM^) - fdx^j , 

where ipj are Hecke eigenfunctions and iV — > oo through primes. For 
this purpose they introduced the quadratic form Q(n) = u(nA, n) (with 
uj(n, m) = n\mi — n^vnx the standard symplectic form) and used it to 
define twisted Fourier coefficients. For a smooth function / e C°°(T 2 ) 
with Fourier coefficients f(n) for n 6 Z 2 , the twisted coefficients are 
given by 

(1.2) f*{u) = (-i) nin2 /w- 

Q{n)=v 

Conjecture (Kurlberg- Rudnick [IS])- As N — > oo through primes, the 
limiting distribution of the normalized matrix elements Fj N ^ is that of 
the random variable 
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where U u are independently chosen random matrices in SU(2) endowed 
with Haar probability measure. 

As evidence, the second and fourth moment were computed to show 
agreement with this conjecture. In particular, the moment calcula- 
tion implies that the limiting distribution is not Gaussian, in contrast 
to generic chaotic systems where the fluctuations are believed to be 
Gaussian [5J [7] . 

In this paper we further study the matrix elements for the cat map 
for composite N. In fact, it is sufficient to understand the case of prime 
powers (see [T^l Section 4.1]), and so we restrict ourselves to this case. 
For N a prime power, we give an explicit formula for the matrix el- 
ements as a weighted sum of certain exponential sums. We then use 
this formula to show that there are sequences of eigenf unctions such 
that the matrix elements decay like iV -1 / 3 rather then the expected 
rate of N~ l l 2+e . We further show that when N = p k with k > 1, the 
matrix elements have a limiting distribution as p — > oo. This distri- 
bution is not Gaussian and it is also different from the (conjectured) 
distribution for k = 1. Instead of behaving like traces of random ele- 
ments from SU(2), here the normalized matrix elements vanish for half 
of the eigenfunctions and for the rest they behave like 2 cos(#) where 
the angle is chosen at random. 

1.1. Results. For every N = p k denote by 

C(p k ) = {Be SL(2,Z/p fc Z)|AB = BA (mod p k )} , 

the group of Hecke operators. For v G Z and x a character of C(p k ) 
define the exponential sum 

E p k(u,x)= e P k ( ux )x(P(x)), 

where 

X(p k ) = {xe Z/p fc Z|(Tr(A) 2 - 4)x 2 ^ 1 (mod p)} 

and 13: X(p k ) "—>■ C(p k ) is an injection of X(p k ) into C(p k ) given by a 
rational function (defined by (13. ip ). 

Theorem 1. For each prime power p k , there is a subset Co(p k ) C C(p k ) 
of characters, with lim^oo ^ Co ^ - = 1 such that 

(1) For any x £ C (p k ) there is a unique Hecke eigenfunction ip, 
s.t., x is a joint eigenvalue. 



4 



DUBI KELMER 



(2) For this eigenf unction, and any elementary observable f n (x) = 
exp(27rm • x) with Q(n) ^ (mod p) 



where Xo is a fixed character of C(p k ) and the sign (±) depends 
on p, k but not on ip. 

If we consider nontrivial prime powers (i.e., k > 1) we can use el- 
ementary methods to evaluate these sums. In particular we find that 
there are matrix elements that decay much slower then the expected 
rate of N~z +e , 

Theorem 2. There are smooth observables f G C°°(T 2 ) ; and sequences 
of Hecke eigenf unctions satisfying \ (0~p N ( f)ipj,ipj) — f r2 f\ N~s. 

We note, however, that these exceptional matrix elements are quite 
rare, in the sense that for a fixed observable the number of matrix ele- 
ments decaying slower then N~2 +e is bounded by 0(p ) (see Corol- 
laryd]) . 

Remark 1.2. In [IT] Olofsson studied the supremum norm of Hecke 
eigenfunctions for the quantized cat map. He showed that for composite 
TV the supremum norm can be of order N*, which is much larger then 
the case of N prime (or square free) where all Hecke eigenfunctions 
satisfy Hoc ^ N e [Him]- Although the two phenomena look similar, 
there does not seem to be any apparent connection between them. At 
least in the sense that the eigenfunctions with large matrix elements 
are usually not the eigenfunctions with large supremum norm. 

For nontrivial prime powers, we can also show that the exponential 
sums E p k{v,x) (and hence also the matrix elements) have a limiting 
distribution as p — > oo. (See [13] for similar results on twisted Kloost- 
erman sums). To simplify the discussion we will assume from here on 

that the observable / is a trigonometric polynomial and let Fj be 
the normalized matrix element as in (11.11) . Let \i denote the measure 
on [0, 71") defined by 



Theorem 3. Let f be a trigonometric polynomial. For any k > 1, 
as p — > oo through primes, the limiting distribution of the normalized 
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matrix elements F) p is that of the random variable 

Y f = 2j2f*(") cos(0„) 

where 9 V are independently chosen from [0, it) with respect to the mea- 
sure /i. 

Remark 1.3. As mentioned above, there can be exceptionally large 
matrix elements for which Fj N ^ ^> N 1 ^ are not bounded. Such matrix 
elements would cause the moments (above the 6'th moment) to blow 
up as iV — > oo. Nevertheless, since the number of exceptional matrix 
elements is of limiting density zero, they do not influence the limiting 
distribution (see section [231 for more details). 

1.2. Outline. The outline of the paper is as follows: In section [2] we 
provide some background on the cat map and its quantization and 
on the notion of a limit distributions. In section [3] we compute the 
formulas for the matrix elements proving Theorem [TJ In section S] 
we compute the exponential sums appearing in these formulas for non 
trivial prime powers, and establish the limiting distribution as p — > oo. 
Then in section [5] we deduce both of the results on the matrix elements 
(Theorems [2] and [3]) from the analysis of the exponential sums. 



2. Background 

The full details for the cat map and it's quantization can be found 
in [15]. We briefly review the setup and go over our notation. 

2.1. Classical dynamics. The classical dynamics are given by the 
iteration of a hyperbolic linear map A G SL(2, Z). 

J G T 2 h-> Ax (mod 1). 

Given an observable / G C°°(T 2 ), the classical evolution defined by A 
is / !-»• / o A. 

2.2. Quantum kinematics. For doing quantum mechanics on the 
torus, one takes Planck's constant to be 1/N, as the Hilbert space 
of states one takes TLn = L 2 (Z/NZ), where the inner product is given 
by: 

yeZ/NZ 
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For n = (711,712) G Z 2 define elementary operators Tjv(n) acting on 
ip G Hn via: 

(2.1) T N (n)ip(y) = e2N(n 1 n 2 )e N (n 2 y)ij(y + ni), 

where e^ix) = e^ 2 ^ . For any smooth classical observable / G C°°(T 2 ) 
with Fourier expansion f(x) = J2 n ^z 2 f(n)e 2mn ' x , its quantization is 
given by 

Op N (f) = £ f(n)T N (n). 

nez 2 

2.3. Quantum dynamics: For any A G SL(2,Z), we assign unitary 
operators Un(A), acting on L 2 (Z/A/"Z) having the following important 
properties: 

• "Exact Egorov": For A = I (mod 2), and any / G C°°(T 2 ) 

^(A)- 1 Pjv (/)^(A) = Op N (f o A). 

• The map A i— » [^(A) is a representation of SL(2, Z/ATZ): If 
C = (mod AT) then U N (A)U N (B) = U N (C). 

We will make use of the following formula for Un(A), (valid for odd 
N and any A G SL(2,Z)) [El Proposition 1.4]. 



where ctat(v4) = Tr(t/jv(-4)) is the character of the representation, 

| ker N (A -J)|=#{ne (Z/iVZ) 2 |n(A - J) = (mod N)} , 

and T N {n) = (—l) nin2 T N (n) are twisted elementary operators. 

Remark 2.1. The twisted operators Tjy(n) have the convenient feature 
that 

C^(A)*f JV (n)C^ r (A)=f JV (nA) 
for any A G SL(2,Z/A^Z) (without the parity condition). 

2.4. Hecke eigenfunctions. Let a > a -1 be the eigenvalues of A in a 
(real) quadratic extension K/Q. Then the vectors v± = (c, a ±jL — a) are 
corresponding eigenvectors w±y4 = a v±. Denote by D = Tt(A) 2 — 4 G 
Z + so that y/D = a — a" 1 . Consider the ring D = Z[a] and denote 
by t : O — > Mat (2, Z) the map sending /3 = n + ma i— > i? = n + mi 
(this map is a ring homomorphism as a and A have the same minimal 
polynomial) . 

For any integer N the norm map, Mk/q '■ K* — > Q*, induce a well 
defined map Af N : (Of NO)* -> (Z/iVZ)*. Let C(iV) = kerA^v be 
its kernel, then its image l(C(N)) C SL(2,Z/A r Z) is a commutative 
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subgroup of SL(2, Z/NZ) that commutes with A (mod N). The Hecke 
operators are then {U^{B)\B G t(C(N))}, and Hecke eigenfunctions 
are joint eigenfunctions of U^{A) and all the Hecke operators. 

The eigenvalues corresponding to each Hecke eigenfunction define a 
character \ of C(N) i.e., U^(l((3))i/; = x(P)ip- We can thus decompose 
our Hilbert space into a direct sum of joint eigenspaces Hn = © x "^x> 
parameterized by the characters of C(N). We say that a character 
X appears with multiplicity one in the decomposition when the corre- 
sponding eigenspace is one dimensional. 

2.5. Limit distribution. We recall the notion of a limiting distribu- 
tion for a sequence of points on the line. For each N let {F^}^ =1 
be a set of points on the line. We say that these points have a limit- 
ing distribution Y (where Y : Q — > R is some random variable on a 
probability space Q) if for any segment [a, b] C R the limit 



#{j\a<F^<b} 



lim — — = Prob(F G [a, b}). 

N^oo iV 

From this definition it is strait forward that making an arbitrary change 
in a density zero set of points (i.e., changing Sn points for each N with 
%■—>()), does not affect the limiting distribution. 

An equivalent condition for having a limiting distribution Y, is that 
for any continues bounded function g the average jjYlj di^j^) con ~ 
verges as iV — > oo to J n ^(y(o;))c?a;. Note that the condition that the 
test function g is bounded is necessary unless both the variable Y and 
the points Fj N ^ are uniformly bounded. In particular, if the points 

Fj N ^ are not uniformly bounded then their moments don't necessarily 
converge to the moments of Y. 

2.6. Notation. We use the notation e(x) = e 2mx . For any iVeNwe 
denote by ejv(-) the character of Z/NZ given by e^x) = e(Jr). When 
there is no risk of confusion we will slightly abuse notation and write 
ejv(f ) f° r eAr(a6 _1 ) (where 6 _1 denotes the inverse of b modulo N). For 
example, for odd and a G Z we may write e2jv(a) = (— l) a eiv(|). 



3. Formulas for Matrix Elements 

For N a prime power we give formulas for the matrix elements of 
elementary observables explicitly as exponential sums. When N is 
prime these formulas appeared in [16] (for primes that split in 0) and 
in [12] (for inert primes). 
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We will make use of the following parametrization of the Hecke oper- 
ators. For any integer 1 < / < k we define subgroups C p (k, I) C C{p k ) 
by 

C p (k,l) = {peC(p k ) |/3 = 1 (rnodp*)}- 
For notational convenience we will also define C p (k, k + 1) = {1}. Let 

X(p k ) = {xe Z/p k Z\Dx 2 ^ 1 (mod p)} 

then the map 

/q -i \ of \ + 1 

(3-1) /%) = -7= - 

is a bijection between X(p k ) and C(p k ) \ C p (k,l) with inverse map 
given by x = ^/^f}^ x ^ (mod p fc ) (note that for (3 ^ 1 (mod the 

inverse map is indeed well defined). For every character \ °f C(p k ) 
and any z/ G (Z/p fc Z)* we have the exponential sum 

# P *0,X)= e P fc ( z/x )x(/?(a ; ))- 

To prove Theorem [1] we will show that for any n G Z 2 with Q(n) = 
v 7^ (mod p), and for every character \ of C(p k ) that appears with 
multiplicity one, the corresponding matrix element is given by 

(f pk (n)i/;,i/;) = n C( k) Yl e P fc (y)xXo(/^)), 

(where xo is a fixed character of C(p k ) and the sign is —1 when p is 
inert and k is odd and +1 otherwise). We can then take our set Co(p k ) 
to be the set of characters appearing with multiplicity one. This set is 
of order p k if p is inert (Lemma 13. 2j) and of order p k — p k ^ 1 if p splits 



Co( P k ) 

£ 

p K 

matrix elements separately for the inert and split cases. 



(Lemma 13. ip . Hence, indeed Co ^ = 1 + 0(^). We will compute the 



3.1. Split case. When p is split, we can give explicit formulas for the 
Hecke eigenfunctions and use them to compute the matrix elements. 
Since we assume that p splits in D, there is a matrix M G SL(2, r Ljp kr £) 

satisfying that M~ l AM = y _i (mod p ). Consequently, the 

Hecke group is given by 

C(p k ) = ^M(^ M- l \x G (Z//Z)*}, 
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where Xo is a fixed character of (Z/p fc Z)* [15J Section 4.3]. 

Lemma 3.1. For any character x o/(Z/p fc Z)* (extended to a function 



normalized joint eigenfunction of all Hecke eigenf unctions with eigen- 
value xXo- Furthermore, if x is not trivial on the subgroup C p (k, k — 1) 
then this is the only eigenfunction. 

Proof. The first assertion is an immediate consequence of (13.21) . For the 
second part, assume that ip is an eigenfunction with eigenvalue xXoi 
and that x is n °t trivial on C p (k, k — 1). Then there is x G C p (k, k — 1) 
with x(xo) ^ 1. Now, let <p = Upk(M)- 1 ^, then for any x G (Z/p fe Z)*, 



hence 4>(xy) = x(x)4>(y). For any y = (mod p) we have that x y = y 
(mod p k ) (as x = 1 (mod^" 1 )). Consequently, (f>(y) = (p(x y) = 
x(xo)4>(y) implying that (p(y) = 0. On the other hand, for y ^ 
(mod p) we have <p(y) = so is uniquely determined (up to 



Remark 3.1. In the case that the character x is trivial on the group 
C p (k, I) (but not on C p (k, I — 1)) then the above argument implies that 
the corresponding eigenspace is of dimension k — I + 1. 

Proof of Theorem [TJ (split case). Let x be a character not trivial 



on C p (k,k — 1). Then, ip = . / -^U p k(M)x is an eigenfunction with 



on Z/p fc Z by setting x(p x ) — ®)> the function vp = */ -^U p k(M)x 




normalization). 



□ 




Consequently, 



for any (twisted) elementary observable 



(r p fc(n)^,^) 



P 



{U pk (M)*T pk (n)U pk (M) X ,x) 



p — 1 



P 



(T p k(m)x,x) 



p — 1 



with m = nM (mod p k 
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Now, let d = y—y 1 so that d 2 = D (mod p k ) (recall Ti{A) = y+y 1 
(mod p k )). Then 



{T pk (m)x,X) = — 2^ e P fc ( 2 ) e P fc ( m 2x)x( ) 

1 r—s dm\m2 dt+1 



p k ^ p \ 2 r 

tex(p fc ) 

where we made the change of variables 2x = rrii(dt— 1). Finally, notice 
that for m = nM (mod we have that 

Q(n) = Lu(nA, n) = mim 2 (t/ — y~ v ) = dm\m2 (mod p k ). 

Hence indeed 

,Q(n)x 



iex(j) fc ) 



□ 



3.2. Inert case. First we show that for p inert, any joint eigenspace 
is one dimensional. 

Lemma 3.2. For N = p k and p inert, the dimension of any joint 
eigenspace satisfies dim7i x < 1. 

Proof. The trace of the quantum propagators satisfy [TJl Corollary 1.6] 

\Ti{U pk {B))\ 2 = ||n6 (Z/p*Z) 2 : n(S - /) = (mod p k )} . 

For p inert, the group C{p k ) is of order #C(p fc ) = p* -1 (p + 1), and the 
groups 

C p {h,i) = {peC{p k )\p = l (modp 1 )}, 

are of order #C p (fc, /) = ^fe) — p fc ~'- Moreover, for any /3 G C p (A;, i) \ 
C p (k,l + 1) we have |Tr(f/ p fc(i( / 9)))| 2 = p 2 '. Consequently 

k-l 

\Tr(U pk m))\ 2 = P k + J2 E |Tr(M^)))lW 

/3eC( P fc ) J=i /3ec p (k,i)\Cp(k,i+i) 

k-l 

= p k + J2( pk ~ l - p k ~ l ~v + p 2k 
i=i 

= p k +p 2k ~ 1 -p k +p 2k =p k #C{p k ) 
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On the other hand, if we denote by n x = dim Ti x then 



' peep) x 
Comparing the two expressions we get 

£-* = ^ky £ |Tr(^(^)))| 2 =/ = dim^ = ^ 

x /3ec(p fe ) x 

Since n x are non negative integers this implies n x < 1. □ 

After establishing this fact, following the idea of Gurevich and Hadani 
[9], we can write the matrix elements of elementary observables as 

(Tpkin)^,^) = Tr(T pk (n)V Xj ), 

with 

the projection operator to the (one dimensional) eigenspace spanned 
by ipj. We then use formula (12.21) for U p k(i(f3)) in order compute 
Tr(T p fc(n)[7pfc(i(/3))). However, in order to do this we first need to give 
a formula for the character of the representation a(B) = Ti(U p k(B)) 
(appearing in (12.21 )). for any B G i(C(p k )). 

Proposition 3.3. There is a character xo £ C(p k ) such that for any 
(3 G C{p k ), we have 

Tr(U pk m)) = (-l) k (-p) l Xo(f3), 
with 1 < I < k the maximal integer such that (3 = 1 (mod p l ). 
Proof. For any 1 < I < k + 1 consider the subgroup of characters 
C<V) = \ x e C(p k )\ x ((3) = 1, V/3 e C p (k,l) 



For 1 < I < k, the group C^ l \p k ) is the kernel of the restriction map 
from C(p k ) to C p (k,l) and hence of order ^qh, ^ = p l ~ l {p + 1) (and 

for / = k + 1 we have C i - k+1 \p k ) = C^{p k ) = C{p k )). 

We will first prove the following: For each 1 < I < k + 1 there is a 
character xi £ C^ l \p k ) and a subset Si C C*^ _1 )(p fc ) of order #5; = p l ~ x 
such that for any (3 G C(p fc ) \ C p (k, /), 

Tr(CV W))) = (~l) k+l+1 XiXi + i ■ ■ ■ XM 

X&Si 
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First for I = k+1 we take the character to be the trivial character and 
the set S k+ i C & k \p k ) = C(p k ) to be the set of characters that appear 
in the decomposition of Ti, p k (there are p k such characters each appear- 
ing with multiplicity one). Then indeed Tr(C/ p fc(i(/3))) = J2 xeSk+l x(P)- 

If k — 1 the sum is over all but one of the characters, say Xo £ C(p), 
and hence Tr (U p (l(/3))) = —Xo(P) as claimed. For k > 1 we proceed 
by induction as follows. 

We assume the assertion is true for 1 < I < k + 1 and show that it 
is true for I — 1. For any f3 G C(p k ) \ C p (k, /), by our assumption 

Tr(U pk (i(P))) = (-l) k+l+1 XiXi + i---xMj2*(P) 

xes ; 

with Si C C {l - 1) (p k ) of order #5, = p'" 1 . The order 4C^{p k ) = 
p 1 ^ 1 + p l ~ 2 hence the complement Sf in C^ _1 )(p fe ) is of order p l ~ 2 . 
Now, if P G" C p (k, I — 1) then the sum over all characters in C ( ^ l ~ 1 \p k ) 
vanish, and hence Z) x es, x(P) — ~ J2 x es c x(P)- We thus have that 

Tr(U pk (l((3))) = (-l) k+l XiXi + i---xM 

xeSf 

On the other hand, for /3 G C p (k, I — 2) \ C p (k, I — 1) we have that 
|Tr(C/pfc(i(/3)))| = p l ~ 2 , which could happen only if x(P) takes the same 
value for all x £ S'f- Now take xz-i to be any character from Sf and 
let Si-! = xT-iSf- Then C C^tf) is of order p'" 2 and 

Tr(f/ pfc (^))) = (-l)^-ViX^ + i---X fc (/3) £ 

Now, let Xo — Xi ' X2 ■ ■ ■ Xk and let f3 G C P (A;, I) \ C p (k, I + 1). Since 
(3 G" C p (/c, / + 1) we have, 

Tr(U pk (L(m = (-l) k+l Xi + iXi + 2---Xk(P) £ 

XG5 ! + i 

On the other hand we also assume /3 G C p (k,l), hence, for all x £ 
Sj+i C (7 w (p fc ) we have x(P) = 1 implying that J2 x es l+1 x{P) = 
#S l+1 = p'. Also for any m < Z, e C^G^) C CW(p fc ), so 
Xm(P) — 1- We thus get that indeed 

Tr(f/ pfc (^))) = (-l) fc (-p)^ (/9). 



□ 
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Proposition 3.4. Let n G Z 2 and B G i(C(p k )). For B = I (mod p) 
the trace Tr{T pk {n)U pk (B)) = 0. Otherwise, there is x G X{p k ) such 
that B = l((3(x)) and 

Tr{f pk {n)U pk {B)) = (-l)^ (/3(x))e p ,(-^^). 
Proof. Use formula (12. 2p for U pk (B) to get that 
Tr(f pk (n)U p k{B)) = 

\kerAB B - I)\p* £ MfAn)f pk ( m )f pk (-rnB)) 

Note that up to a phase T pk (n)T pk (m)T pk (— mB) = e ia T pk (n—m(B—I)) 
and recall that Tr(T pk (n)) = unless n = (mod (see e.g., [TB] 
Lemma 4]). Hence, the only summand that does not vanish is the one 
satisfying n = m(B — J) (mod We can assume n^O (mod p), so 
that the trace vanishes whenever B = I (mod p). Otherwise, B = t{j3) 
for some (3 G C(p k ) \ C p (k, 1) and cr pk (B) = (— l) fc xo(/3) so that 

Tr(f pk (n)U pk (B)) = (-l) fc Xo (/?)e pfc (-^^), 

with m = n(B — J) -1 (mod 

Now recall the parametrization C(p k )\C p (k, 1) = |/3(a;) |x G X(p fe )}, 

with (3(x) = ^) X x + \ - We claim that for B = l((3(x)) and m = nB we 
have that u(m,mB) = Q{n)x. To show this substitute ((3(x) — = 
and (p(x) - l) _1 /3(a;) = ^%±±. Consequently we get 

x 



uj(m, mB) = uj ( ),m( ^ ) J = — Lu(m(V D),n). 

Recall that \f~D = (a — o _1 ) so that indeed 

u}{ni{s/D), n) = uo(n(A - A' 1 ), n) = 2lu(uA, n) = 2Q(n). 

□ 

Proof of Theorem [TJ (inert case). For every character x let 
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be the projection operator to the (one dimensional) eigenspace corre- 
sponding to x- Let ip be the corresponding Hecke eigenf unction. Then 

CZ>(n)V,V> = Tr(T>(n)P x ) = — ^ £ Tr(T pfc (n)C/ pfe (S))x(S). 

j fleC(p») 

Now from the Proposition 13.41 

Tr(f p k(ri)U p k(B)) = (-lf Xo (B)e pk (-^^), 
implying that 

After a change of variables x i— ► — x we get 

(-l) fc 

(T p k(n)ij,ip) = ^ c Epk(Q(n)/2,xXo) 



□ 



4. Analysis of the Exponential Sums 

In this section we compute the exponential sums E p k(v, x) for any 
prime power k > 1. This can be done using elementary methods (see, 
e.g., [11, section 12.3] or [TJ Chapter 1.6]), however, since the setup 
here is slightly different we will perform this computation in full. We 
then evaluate all mixed moments of these exponential sums to deduce 
their limiting distribution. 

4.1. Computation of exponential sums. For v e Z/p k Z its "square 
root" (modulo p k ) is the set 

Sq(v } p k ) = {xe Z/p k Z\x 2 = v (mod p k )} . 

Note that for v ^ (mod p) this set contains two or zero elements, 
for v = (mod p k ) it contains p[ fc / 2 l elements (and for v = p l v with v 
coprime to p it contains zero or 2p l l 2 elements). 

Proposition 4.1. For k = 21 even 

*eSq(^±V) 

Dx 2 ^l(p) 

where t x e Z/p l Z satisfies that x(l + p l y/Dx) = e p i(t x x) 
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For k = 21 + 1 odd 

Dx 2 ^l(p) 

where t x £ Z/p z+1 Z satisfies x(l + p l VDx + p 2l ^x 2 ) = e p i+i{t x x), and 
Q(x) is the Gauss sum given by 

G{x)= Y e p(f( x )v 2 + 9{x)y), 



with f(x) = ^%zr an d 9{x) = p~ l iy — t x D J>_ 1 )- (Notice that for 
x £ Sq( 2t *p U , p l ) we have [y — t X D 2_ 1 ) = (mod /ience — 
t x j5^—i) gives a well defined residue modulo p). 

Proof. First for k = 21, write the sum as 

E p k(v,x)= Y Y e P k ( u ( x + p l y))x(f3(x + p l y))- 



Replace f3(x + p l y) = + ^(x)p l y) (mod p 2i ) to get 





E p k(v,x)= Y e P k ( ux )x(i3(x)) Y e P i {vy)x{ l +p l ^{x)y)- 



xex(p') yei/p l . 



Differentiating (5{x) = we get ^-(x) = — ££P_ X , so that the inner 

sum takes the form 



The map a: i — ^ 1 + p l \fT)x is an isomorphism of Z/p'Z and C p {2l,l). 
Hence, for any character x °f C(p 21 ) there is t x £ TLjp l TL such that 
x(l +p l VDx) = e p i(t x x). We can thus write the inner sum as 

Y e p i{vy)e p i{- ^jf* ^ = E 6 p' ^ ~ ^ _ i ^ ' 



This sum vanishes unless x £ Sq( 2tx +" ; p 1 ') in which case it is equal p l . 



vD 

Now for k — 21 + 1, we start again by writing 



E p k(v,x)= Y Y e P k ( u ( x +p l y))x((3(x+p l y)), 

xGX(p l ) y£Z/p l Z 
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and replace 

P(x + p l y) = P(x)(l + ^{x)p l y + \j(x)p 2l y 2 ) (mod p 2l+1 ). 

It easy to verify that the map x i— > 1 + p l \f~Dx + p 2l ^x 2 is an isomor- 
phism of Z/p l+1 Z with C p (2l + 1,1). Consequently, for every character 
X of C(p k ), there is t x G Z/p l+1 Z such that 

x(l+p l VDx+p 2l jX 2 ) = e pl+ i(t x x). 

By differentiating (5{x) = ^^+* (twice), we get that 

= ^ { _2j^^ )+pl D { 2j^ 1)2 (modp , +1) _ 

implying that the inner sum is of the form 

2t x 

Dx 2 - 1 V 4 * Dx 2 - 1 



This sum vanishes unless x G Sq( 2 ^ i/ ,p ; ). To see this make a change 
of summation variable y i— > y + p to get that 

yei/p l+1 i 

Now unless z/ — = (mod p') we have that e p i{{v — j^rti )) ^ 1, 

implying that the sum must vanish. For a; G Sq(^=^,p') the inner 
sum given by p l times the Gauss sum 

G( x )= e p(f( x )v 2 + 9{x)y)- 

□ 

In particular this computation implies that for most characters the 
exponential sum has square root cancelation. 

Corollary 1. For any character x with 2t x ^ —v (mod p) there is 
9 = 9(x, v) G [0, 7r) such that E pk (v, x) = P k/2 cos(9(v, x))- 
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Proof. The condition 2t x ^ —v (mod p) implies 2t *p U 7^ (modp). 
Hence for any 1 < I < k 

nc ( 2t x + v h_ / 2 2 -^ = U (modp) 
* q[ vD ,P) ~ \ " otherwise 

Now for k = 21 even, recall that 

xesq(^y) 

If ™xgL ^ □ (mod p) this sum vanishes. Otherwise it is a sum over 

two elements of absolute value p l = p k l 2 hence indeed E pk (v,x) = 

2p k / 2 cos(9(u, X )). 

For k — 21 + 1 odd we have 

E p k(is,x)=P l e pk (vx)x((3(x))G(x). 

The condition 2t x ^ — i/ (mod p) implies that the Gauss sum G(x) is 
not a trivial sum and hence of order y/p. As before, the sum E p k(u, x) 
either vanishes (if 2t *p V ^ □ (mod p)) or it is a sum of two elements 
of absolute value p l+ ^ = p k l 2 . □ 

On the other hand, if 2t x = — v (mod p 21 ') for some V < | then 
the sum contains p z ' elements and could be much larger. Moreover, in 
the odd case, if 2t x = —v (mod p l+1 ) then g(x) = (mod p). Also, 
in this case any x G Sq( ^jy- ,p l ) satisfies x = (mod p) and hence 
also f(x) = y^fzi = (mod p). So that in this case the Gauss sum 
\G(x)\ = p rather then ^Jp and the sum is even bigger. In particular 
we get 

Corollary 2. For and x £ ^(p 3 ) with t x = — z/ (mod p 2 ) 7 

l^(^,x)l =p 2 - 

4.2. Equidistribution of exponential sums. We now show that 
as p — > oo the normalized exponential sums p~ k ^ 2 E p k{y ) x) become 
equidistributed with respect to the measure 

(4.i) M ) = \ f fy + h f fWe - 

— — J Q 

For fixed v and a character x-, if ^T£r~ i s n °t a square modulo p then 
the sum E pk {v,x) — (or equivalently 0(v,x) — f)- The following 
lemma shows that this happens for roughly half the characters, and 
that this behavior is independent for different values of v. 
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Lemma 4.2. Fix a finite set v = {ui, . . . , u r } of nonzero distinct inte- 
gers. Then, 

-# It G Z/pZlVj, ^ = □ (modp)) = ^ + 0(^). 

Proof. We can write 
2 r # |t g Z/pZ|Vj, = □ (mod p) 



En(»( 



t j=i 

with X2 the quadratic character modulo p. Now expand the right hand 
side 

* j=l v 3 7 JC{i,...,r} * ieJ 3 

Where the sum is over all subsets J C {1, . . . , r}. The contribution 
of the empty set is exactly 1 = p, while for nonempty J we get an 
exponential sum of the form J^Xadljej I^T 1 )- Since we assumed all 
Uj are distinct, the polynomial g(t) = Yljej i^r * s no ^ a sc L uare an d we 
can apply the Weil bounds J2tX2(g(t)) = 0( v /p) [18J. Consequently, 
we have that indeed 

2 r # jt G Z/pZ|Vj, ^^ = D (modp)J =p + 0(Vp). 

□ 

Next we need to show that for the rest of the characters (when the 
exponential sum does not vanish) the angles 9(x, v) become equidis- 
tributed (independently) in [0, 7r]. We will do that by computing all 
mixed moments. However, we recall that there are exceptional charac- 
ters for which the normalized exponential sums are not bounded caus- 
ing the moments to blow up. For that reason we first restrict ourself to 
a set of "good" characters (of limiting density one) for which the sums 
are bounded and only then we calculate the moments. 

Fix a finite set of r nonzero distinct integers v = {z/ 1; . . . , z/ r }, and 
define the set of "good" characters to be 

S+i?) = {xe C(p k )\Vj,2t x ± - Vj (mod p)} , 

where t x is determined by \ as above. Then for any character x £ 
S p k(u), we can write E p k(uj,x) = p fc//2 cos(6(vj, x)) with 0(vj,x) e 
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[0, 7t). Furthermore, for any Vj there are precisely p k ~ 2 (p±l) characters 
with 2t x = Vj (mod p) (this is the size of the kernel of the restriction 

map from C(p h ) to C p (k, 1)). Hence, ^ p ^ -I = 1 + O(^) and the set 

S p k(v) is of (limiting) density one inside C{p k ). 

Before we proceed to calculate the moments we will need to set some 
notations. For any k define the set 

Y(p\ v) = {x G (Xip^flu^Dxl - 1) = Vj{Dx) - 1), V2 < j < r) 

For every fixed set of integers n = {ni, . . . , n r } let 

Y (p k ,v,n) = jz G Y(p k )\ Y[P(xj) n * = 1 (mod p k 

For notational convenience we will sometimes use the notation S p k, 
Y(p k ), Y Q (p k ) where the dependence on v and n is implicit. We will 
also denote by Y'(p k ) (respectively y o '(p fc )) the elements of Y(p k ) (re- 
spectively Y (p k )) with all Xj ^ (mod p). 

Lemma 4.3. As p — > oo, the number of points in Y'(p k ) satisfy 

#Y'{p k )=p k + 0(p k - 1 i) 

Proof. For any t G TLjp kr L satisfying Vj, t ^ Uj (mod p) we have that 

On the other hand if t = Vj (mod p) for some j, then 

#{x|Vj, v J {Dx)-l)=t}<2 r p k - 1 

(as there are at most two possibilities for x« with i ^ j and at most 
2p k ~ 1 possibilities for Xj). We thus have 

= 2y- x # ji G (Z/pZ)* |Vj, = □ (mod p) J + o^- 1 ). 

Also note that #Y'(p k ) = #Y(p k ) + 0(p k ' 1 ). To conclude the proof 
we use the estimate 

2 r # jt G Z/pZ|Vj, = n (modp)J=p + 0(v^), 

from lemma l4~2l 

□ 
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Lemma 4.4. As p — > oo, the number of points in 5^(p fc ) satisfy 

Wo(p k ) = oip 1 *- 1 ). 

Proof. To prove this bound we will show that there is a nonzero poly- 
nomial F(t) with integer coefficients such that for any x G Y (p k ), with 
b = ^li£^~ll we have = (mod p k ) (recall that for x-y G X(p k ) 

we have Dx\ ^ 1 (mod p) and hence b ^ (mod p) is invertible). This 
would imply that b can take at most deg F values modulo p, implying 
that Wo{p k ) < 2 r deg{F)p k -\ 

Now to define F, consider the formal polynomial in the variables 
Pf 1 ,- ■ -fit 1 given by 

G{^...(3 r ) = J] m j /" 1 ) • 

<re{±i}'- \i=i / 

Recall that if a polynomial in two variables x, y is symmetric under 
permutation then it can be written as a polynomial in the symmetric 
polynomials o\ = x + y, a 2 = xy (see e.g., [3[ Chapter 6]). The poly- 
nomial G is symmetric under any substitution (3j i— > (3J 1 and hence 

there is another polynomial F in r variables with integer coefficients, 
satisfying 

G(f3 1 ,...,f3 r )=F(f3 1 + fc\...,p r + f3; 1 ). 

Define the polynomial F(t) = F(2 + uxt, . . . , 2 + u r t). For any 
Xx,...,x r with x 2 j = 2h y p (mod p k ) we have ^(a^-)-!-^^) -1 = 2+Vjb~ 1 

(modp k ) (recall f3(x) = ^^zj)- Hence, 

G^an), . . . , P(x r )) = F(2 + v x b~\ ...,2 + u r b~ l ) = Fib' 1 ). 

Now, if in addition (3(xi) ni ■ ■ ■ (3(x r ) nr = 1 (mod p k O) then indeed 
F^- 1 ) = G(J3(xi), /3(x r )) = (mod p k ). 

It remains to show that F(t) is not the zero polynomial. To do this, 
we think of it as a complex valued polynomial, and note that for it to 
be identically zero there has to be some choice of signs a G {±l} r so 
that the function 
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satisfies G a {t) = 1. Assume that there is such a choice a, so the 
derivative G' a (t) must also vanish. But we have 

r 



so as t — ► — y the term t 2^ t ]_ u ^ blows up while the rest of the terms 

remain bounded (recall that all Vj are different). In particular G' a (t) is 
not identically zero. □ 

Remark 4.1. The bound ^Yq{p) = 0{p k ~ l ) is probably not optimal. 
Notice that if the polynomial F(t) defined above is separable (i.e., if it 
has no multiple roots) then there are at most deg F solutions to F(t) = 
(mod p k ) and the corresponding bound would be i^Y^ij^) = 0(1). 

We now preform the moment calculation establishing the limiting 
distribution of the exponential sums (when running over characters in 
S p k). 

Proposition 4.5. Let \i be as in &4-l\) and let g G C([— 1, l] r ) be any 
continuous function then 



lim ~k 52 9(cos(9(ui,x),---,cos(e(u r ,x)) = 

s pk 

g(cos(9i), . . . , cos(0 r ))d/j,(9i) ■ ■ ■ dji{9 r 



Proof. We will give the proof for k = 21 even, the odd case is analogous. 
Since we can always approximate the function g by polynomials, it is 
sufficient to show this holds for all monomials of the form 

g{x) = (2x 1 ) m ^--(2x r ) m r. 

We thus need to show that 

p^oo p k ^ V ./ro-1 



With out loss of generality we can also assume that all the rrij are 
nonzero (since /i is a probability measure, if rrij = then the corre- 
sponding factor is 1 and we can consider the same problem for r — 1 
instead of r). In this case the right hand side is given by 

J! QT (2cos(0)r d^(9)) = II Q [ (2cos(9)r f ) • 
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The integral in each factor is ~ ( m ^ | for m~ = 2n,- even and it is zero 

2 \n 3 ) i 3 

otherwise. 

Now fix a character \ £ and let t x £ Z/p'Z as above. If 
Sq( 2 ^^ ,^) = then 2cos9(uj,x) = p~ k / 2 E p k(b'j,x) = 0. Otherwise, 

2cos(6(v vX ))=p- k/2 E pk ^,x) = 2M(e pk ( 1 ^)x{P{x j ))), 

with £ Sq(^jp-,p l ) (recall that for x £ S p k we know 2t x + v ^ 
(mod p)). Hence, the only contributions to the sum 

n( 2cos ^>x)r, 

comes from characters x such that for all j there is Xj £ Sq( 2< *^ J , p') 
(equivalently, there is a;.,- £ (Z/p z Z)* satisfying Uj(Dx 2 — 1) = 2t x 
(modp')). Also note that if we multiply x by an Y character that is 
trivial on C p (k,l) this does not change t x . Let C^{p k ) be the group 
of characters that are trivial on C p (k,l), and for any b £ Z/p'Z let 
%6 £ C(p k ) be a representative of C(p k )/C^(p k ) with t X6 = 6. We thus 
have that 

l £ n(2cos^,, X )r = 
= 2^ £ £ n( 2 «» ^.xxor 

iey(p ; )xec , (')( p fe) j 

where 6 = 6(x) = ^(^i x ) . 
Now use the formula, 

(2cos(0)r = f;H cos((ro-2n)0). 

n=0 ^ ' 

The main contribution comes from the terms where in each factor rrij — 
2rij = 0. This vanishes unless all rrij are even in which case it is given 
by 

Lee nfc)=nifcV^ 



where we used Lemma 14.31 to get that #Y'(p l ) ■ ^C^(p k ) = p k + 
0(p fc ~5). 
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It thus remains to bound the rest of the terms, which is reduced to 
the vanishing (in the limit p — > oo) of the sums 

for any nonzero integers {n\, . . . ,n r }. 
For any x G Y'(p l ) we have that 



cos(V(^XX*) = 2R{e pk {^p-)xxm^) ni )) 

with b = 8/1 x ) _ When expanding the product Ylj cos(rij9(vj, xXb)) 
we get a sum over 2 r terms, each of the form 

3=1 3=1 

We thus need to bound the exponential sum coming from each term. 
We will now bound the corresponding sum 



i E E v(E^)w 4 (f[-3fe)" j 

5;ey'(p i )xeC' (i) (p fc ) ■ 7=1 ■ ?=1 



(the same bound obviously holds when changing any rij to — rij). Now, 
rewrite this sum as 

ier'(p') 3=1 3=1 xec(')( p fc) 3=1 

and note that the inner sum vanishes unless Ylj=i 0{ x j) nj = 1 (mod p') 
in which case it is equal ^C^(p k ) = p k ~ l . We can thus rewrite this 
sum as 

S&^p 1 ) 3=1 3=1 

which is trivially bounded by p~'#E '(p') = O(-) (Lemma 14.41) . □ 

Remark 4.2. The above proof also gives the rate at which the fluc- 
tuations of the normalized exponential sums approach their limiting 
distribution. If one takes the test function g in Proposition 14.51 to be 
smooth then the rate of convergence is 0(— p). This rate comes from 
the bound on the error term in Lemma 14.31 which seems to be a sharp 
bound. 
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5. Back to Matrix Elements 

We can now deduce Theorems [2] and [3] from Theorem [I] and the 
analysis of the exponential sums. 

Proof of Theorem [2]. Let f(x) = e 2mn ' x be any elementary observ- 
able. Take N = p 3 to be a prime cubed. Then by Corollary [2] there is a 
character satisfying \E^(^^-, x)\ — P 2 - Let ^ be a Hecke eigenfunction 
corresponding to \i then by Theorem [T] we get 

□ 

Proof of Theorem [31 Let / be a trigonometric polynomial and write 

\n\<R 

for some fixed R > 0. Let {vx, . . . , is r } = {Q(n)\0 < \n\ < R}, and 
consider the random variable 

r 

Y) = 2^/#(^)cos(^). 

3=1 

with 9j chosen independently from [0, 7r) with respect to fi. We need 
to show that as p — > 00 the limiting distribution of F^ p ■* is that of Yf. 

For any character % of C(p k ) consider the weighted sum of the cor- 
responding exponential sums 

By Proposition 14.51 as p — > 00 the limiting distribution of i 7 ^ as x 
runs through 5 p fc (hence, also as x runs through the whole group of 
characters) is that of Yf. 

Now, for p sufficiently large (i.e., p > max{z/j}) and Xj £ C (p k ), we 
have 

If we further assume that XjXo £ S p k(u) then |£ , p fc(y, xXo)| < 2p fc//2 , 
and hence 

ff> = ^f> + (i). 
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The set of characters {xj £ Co{p k )\xjXo £ S p k} is again of density 

one, hence, the limiting distribution of F- is the same as of 
concluding the proof. □ 
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